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Abstract 



These two lectures give a pedagogical introduction to the "string- 
inspired" worldline technique for perturbative calculations in quan- 
tum field theory. This includes an overview over the present range of 
its applications. Several examples are calculated in detail, up to the 
three-loop level. The emphasis is on photon scattering in quantum 
electrodynamics. 
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Introduction 

Many one-loop amplitudes in relativistic quantum field theory admit represen- 
tations in terms of path integrals over the space of closed loops [|l], [2], |||, |3], |5|, |6|, • 
Those path integrals are closely related to the Feynman path integrals used in 
nonrelativistic quantum mechanics, though they have perhaps not reached quite 
the same practical importance. Still there exists a considerable amount of lit- 
erature on applications of particle path integrals to quantum field theory (see 
e.g. H H [U], 11, [l3|, 14|), to which I cannot possibly do justice here. 



The subject of these lectures rather is a specific and novel way of evaluating this 
type of path integral. This method may be called "string-inspired" because it is 
analogous to calculations in string perturbation theory, and its development was 
triggered by efforts [[H], [l(| to use the peculiar organization of string amplitudes 
to improve on the efficiency of calculations in ordinary quantum field theory. 
However this is to some extent accidental, and the practical application of this 
technique requires no knowledge of string theory. 

The lectures are meant as an introduction for practicioners in quantum field 
theory. More effort will therefore be spent on explaining the calculation rules 
of the worldline formalism, than on rigorous derivations. 

The first lecture is devoted to the theoretical side. I will first sketch the history 
of the subject, and then give an overview over the range of one- loop amplitudes 
in field theory which are presently accessible to this method. I concentrate 
then on the explanation of an even more recent multiloop generalization of this 
formalism ||§ ||| g g| H H> @ • 



The second lecture contains four detailed examples of calculations performed 
with this formalism, mainly taken from quantum electrodynamics: 

1. One - loop QED vacuum polarization. 

2. One - loop QED photon splitting. 

3. The two - loop (scalar and spinor) QED (5 - functions. 

4. The three - loop scalar master integral. 



1. First Lecture: Theory 

1.1) The Bern-Kosower Formalism 

It has been realized not too long ago that string theory may be used as a 
guiding principle for finding new and perhaps more efficient parameter integral 
representations for amplitudes in ordinary quantum field theory. The relevance 
of string theory for this purpose resides in the well-known fact that, in a sense 
which can be made precise, string theory contains ordinary field theory in the 
limit where the tension of the string becomes infinite. In this limit, all massive 
modes of the string get suppressed, and one remains with the massless modes. 
Those can be identified with ordinary massless particles such as gauge bosons, 
gravitons, or massless spin-^ fermions. 

To use this fact for actual calculations of amplitudes may appear an enormous 
detour. It is worthwhile nevertheless due to the different organization of string 
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amplitudes, and of the different methods used to calculate them. For details 
on string perturbation theory, the reader may consult |^] ; let us mention here 
just the following qualitative points which will be important later on: 

1. Scattering amplitudes in string theory are usually calculated in first, not 
second quantization. 

2. String theory calculations involve a much smaller number of different 
"Feynman diagrams". 

The second point is illustrated in fig. 1, depicting a two-point "Feynman dia- 
gram" for the closed string. In the infinite string tension limit this Riemann 
surface is squeezed to a graph, though not to a single one - two Feynman 
diagrams of different topologies emerge. 




Figure 1: Infinite string tension limit of a string diagram 
The first actual calculation done along these lines was performed by Metsayev 



and Tseytlin [46] in 1988, who managed to extract the correct one- loop /3 - 
function coefficient for pure Yang- Mills theory from the partition function of an 
open string propagating in a Yang-Mills background. This yielded also some 



explanation for a fact which had been noted long before [47], namely that 



this /3 - function, when calculated in dimensional regularization, vanishes for 
D = 26, the critical dimension of the open string. 

A systematic investigation of the infinite string tension limit was undertaken 



in the following years by Bern and Kosower 15]. This was done again with a 
view on application to QCD, however now to the calculation of complete on- 
shell scattering amplitudes. Again the idea was to calculate, say, gauge boson 
scattering amplitudes in an appropriate (super) string model containing SU(n) 
gauge theory, up to the point where one has obtained a parameter integral 
representation for the amplitude considered. At this stage one performs the 
infinite string tension limit, which should eliminite all contributions due to 
propagating massive modes, and lead to a parameter integral representation 
for the corresponding field theory amplitude. 

For the calculation of the string scattering amplitude, one may use the Polyakov 
path integral. In the simplest case, the closed bosonic string propagating in flat 
spacetime, this is of the form 
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Vi ■ ■ ■ V N ) ~ Yl I Vh I Vx ( a > r ) yi • • • Vn e~ S[x ' h] . (1) 

top •* 

Here the / Vx(a, r) is over the space of all embeddings of the string world- 
sheet with a fixed topology into spacetime. The / Vh is over the space of all 
worldsheet metrics, and the sum over topologies Stop corresponds to the loop 
expansion in field theory (fig. 2). 





Figure 2: The loop expansion in string perturbation theory 
S[x, h] is a free gaussian action, 

S[x,h] = I dadTVhh aP r} lJLV d a x lx d3X v . (2) 

47TCK J 

Here ^7 is the string tension, and the vertex operators V±, . . . , Vn represent 
the scattering string states. In the case of the open string, which is the more 
relevant one for our purpose, the worldsheet has a boundary, and the vertex 
operators are inserted on this boundary. For instance, at the one-loop level the 
worldsheet is just an annulus, and a vertex operator may be integrated along 
either one of the two boundary components. 
The vertex operators most relevant for us are of the form 



V4 = / dre ik - x ^ (3) 



V A 



J dTT a e li i^e ih < T \ (4) 



They represent a scalar and a gauge boson with definite momentum and po- 
larization. T a is a generator of the gauge group in some representation. The 
integration variable r parametrizes the boundary in question. Since the action 
is gaussian, / Vx can be performed by Wick contractions, 



^(riK(r 2 )) = G(r 1 ,r 2 )^, (5) 

where G denotes the Green's function for the Laplacian on the annulus, re- 
stricted to its boundary. 

Performance of the one-loop path integral then leads to the following Bern- 
Kosower Master Formula for the (single-trace partial contribution) to the one- 
loop N-gluon scattering amplitude: 
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JV , fOO jrp N-l 

a N]1 ~ (a')~- 2 z —[AnT]-- n dne(n-T i+1 ) 

Jo T ^ 



x cxp 



Oi'GijPiPj + Va'Gij(p i ej - pj£i) - GijEiEj 

l<i<j<N 



I multi— linear • 

(6) 



Here it is understood that only the terms linear in all the polarization vectors 
£i, . . . ,£jv have to be kept. Dots generally denote a derivative acting on the 
first variable, G(t\,tz) = ^GfYi, T2), and we abbreviate Gij := G(ri,Tj) etc. 
D is the spacetime dimension, and Z denotes the partition function of the free 
string. The T - integral is a remnant of the path integral over the space of all 
metrics. 

This amplitude is still to be multiplied with a colour trace factor tr(T ai • • • T ajv ), 
and to be summed over non-cyclic permutations (in the original Bern-Kosower 
approach colour decomposition is used for the adjoint representation, so that 
only the fundamental representation matrices appear). 

The main point to be noted is that this formula serves as a unifying generating 
functional for all one-loop N-point gluon amplitudes - something for which no 
known analogue exists in standard field theory. 

The analysis of the infinite string tension limit, for which I refer the reader to 



|15, 16 1 , proceeds differently for the cases of a spin-0, spin-^, or spin-1 particle 
circulating in the loop. The result of this analysis can be summarized in the 
"Bern-Kosower Rules", which allow one to construct the final integral repre- 
sentations without referring to string theory any more. The relation of those 
parameter integrals to the ones originating from the corresponding Feynman di- 



agram calculations has been clarified in [17]. Concerning the advantage of the 



Bern-Kosower rules over the Feynman rules, we mention the following points: 

1. Superior organization of gauge invariance. 

2. Absence of loop momentum, which reduces the number of kinematic in- 
variants from the beginning, and allows for a particularly efficient use of 
the spinor helicity method. 

3. The method combines nicely with spacetime supersymmetry. 



This has been demonstrated in both gluon [Ig] and graviton 19, 20] scattering 
calculations. However, we are interested here mainly in a different property of 
this formalism. To demonstrate this property, it suffices to consider the simplest 
possible case of a Bern-Kosower type formula, namely the one-loop N - point 
amplitude for (f> 3 - theory. For this case, which was already considered in ||, 
one has (up to a coupling constant factor) 
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where 



G B {T %1 Tj) =| n - Tj | T 3 . (8) 

This parameter integral can, for any given fixed ordering of the integration 
variables > Tj 2 > ■ • • > Ti N , easily be identified with the corresponding 
Feynman parameter integral [jl5| , |27j (the usual - parameters just correspond 
to the differences Tj fc — Tj fc+1 ). The complete integral therefore does not represent 
any particular Feynman diagram, with a fixed ordering of the external legs, but 
the sum of them (fig. 3): 




Figure 3: Sum of one- loop diagrams in cp 3 - theory. 

This may not seem particularly relevant at the one- loop level. However eq.( |?j) 
holds off-shell, so that one can sew together, say, legs number 1 and N, and 
obtain the following sum of two - loop (N — 2) - point diagrams (fig. 4): 




Now this is obviously interesting, as we have at hand a single integral formula 
for a sum containing diagrams of different topologies. We may think of this 
as a remnant of the "less fragmented" nature of string perturbation theory 
mentioned before (fig. 1). 

Moreover, it calls certain well-known cancellations to mind which happen in 
gauge theory due to the fact that the Feynman diagram calculation splits a 
gauge invariant amplitude into gauge non-invariant pieces. For instance, to 
obtain the 3- loop j3 - function coefficient for quenched QED, one needs to 
calculate the sum of diagrams shown in fig. 5. 

Performing this calculation in, say, dimensional regularization, one finds that 
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+ 




+ 



+ 



Figure 5: Sum of diagrams contributing to the 3-loop QED (3 - function. 



1. All poles higher than - cancel. 

2. Individual diagrams give contributions to the (3 - function proportional 

to £(3) which cancel in the sum. 



The first property is known to hold to all orders of perturbation theory [48|. 
Recently concepts from knot theory have been used to establish a link between 
both properties [49|, and thus to predict the rationality of the quenched (3 - 
function to all orders. 

It seems therefore very natural to apply the Bern-Kosower formalism to this 
type of calculation. However, in its original version the Bern-Kosower formalism 
was confined to tree-level and one-loop amplitudes. A multiloop generalization 
along the original lines has not yet been constructed, though partial results exist 



1 21, p2| . In the absence of better ideas one could, of course, always start at the 
one-loop level, and pursue the explicit sewing procedure indicated above. This 
turns out to be unnecessarily clumsy, though; there is a more efficient way of 
inserting propagators, which will be explained in the third part of this lecture. It 
is based on a purely field-theoretical approach to the Bern-Kosower formalism, 
which was, at the one-loop level, proposed by Strassler [p3fl . This approach uses 
a representation of one-loop amplitudes in terms of first-quantized worldline 
path integrals, of a well-known type. I will first demonstrate this method for 
the example of photon scattering in scalar QED. 

The one-loop effective action for a Maxwell background induced by a scalar loop 
may be written in terms of the following first-quantized particle path integral 
(see e.g. g [H): 



T[A] 



dT 
~T 



-m 2 T 



J Vxexp 



T 

dr 



1 



(9) 



In this formula T is the usual Schwinger proper-time parameter, and m the mass 
of the particle circulating in the loop. At fixed proper-time T, we have a path 
integral over the space of periodic functions x^{t) with period T, describing all 
possible embeddings of the circle with circumference T into spacetime. We will 
generally use dimensional regularization, and therefore work in D dimensions 
from the beginning. The spacetime metric is taken Euclidean. 
The essence of the method is simple: We will evaluate this path integral in a one- 
dimensional perturbation theory. If we expand the "interaction exponential" , 
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cxp 



dr ieA^x* 1 



oo t ■ \N N 



(10) 



Af=0 " i=l ■ 

the individual terms correspond to Feynman diagrams describing a fixed number 
of interactions of the scalar loop with the external field (fig. 6): 




+ 




+ ... 



Figure 6: Expanding the path integral in powers of the background field 

The corresponding N - photon scattering amplitude is then obtained by spe- 
cializing to a background consisting of a sum of plane waves with definite po- 
larizations, 



N 



(11) 



and picking out the term containing every once (this also removes the 
jn in eq.( [U^)). We find thus exactly the same photon vertex operator used 
in string perturbation theory, eqs.(|3|),(||), inserted on a circle instead on the 
boundary of the annulus. But in the infinite string tension limit the annulus 
gets squeezed to a circle; we may therefore think of the path integral eq.(||) 
as the infinite string tension limit of the corresponding Polyakov path integral 
described above. Again the evaluation of the path integral at fixed values of 
the parameters reduces to Wick-contracting 



x ^ e 



X 



HN a ik N -XN 

N e 



(12) 



The Green's function to be used is now simply the one for the second-derivative 
operator, acting on periodic functions. To derive this Green's function, first 
observe that / Vx(t) contains the constant functions, which we must get rid of 
to obtain a well-defined inverse. Let us therefore restrict our integral over the 
space of all loops by fixing the average position Xq of the loop: 



Vx 



x?(t) 



dx jT>y, 
. 



(13) 
(14) 
(15) 
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For effective action calculations this reduces the effective action to the effective 
Lagrangian, r 

T = J d D x C(x ). (16) 

In scattering amplitude calculations, the integral over xq just gives momentum 
conservation, N 

f d D x e ixo ^i k * = (2tt) d 5(jr Jfej) (17) 
J i=l 

The reduced path integral / T>y(r) has an invertible kinetic operator. This 
inverse is easily constructed using the eigenfunctions of the derivative operator 
on the circle with circumference T, {e 2mn T ,n £ Z}, and leads to 

J —2 OO 1 _ rp 

2{T1 |(_) | T2)=2 X^^e-^ = GB(Tl , T2 >--. ( 18 , 

Here Gb is the function which was introduced already in eq.(|8|). The constant 
— T/6 drops out of all calculations once momentum conservation is applied, and 
is therefore usually deleted at the beginning. 

Working out eq.(^) by Wick-contractions using the function Gb leads to a pa- 
rameter integral representation for the one-loop N-photon scattering amplitude 
identical with the one encoded in the Bern-Kosower master formula for this 
special case. We will look at this parameter integral in more detail later on. 
Alternatively, one can also apply this method for calculating the higher deriva- 
tive expansion of the effective action T(A) p4| , p5| , |26fl . This can be done in a 
manifestly gauge invariant way using Fock-Schwinger gauge centered at x$. 

1.2) One-Loop Wordline Path Integrals and Correlators 

The representation of amplitudes in quantum field theory in terms of first- 
quantized particle path integrals is an old and well-studied subject, and I cannot 
possibly survey the literature here. Let me just mention a few contributions 
which were particularly important for the development of the subject: 

1. R.P. Feynman 1950 |Q] (path integral representation of the scalar propa- 
gator). 

2. E.S. Fradkin 1967 (|| (path integral representation of the electron propa- 
gator) . 

3. Berezin and Marinov 1977 | |5Q| , L. Brink, P. Di Vecchia and P. Howe 1977 
[f)l| l (superparticle Lagrangians) . 

4. L. Alvarez-Gaume 1983 |l0| (application of superparticle path integrals 
to anomaly calculations). 

The following overview is restricted to those cases which are presently relevant 
for the "string-inspired" technique. 



S 



The Free Gaussian Path Integral 

With our conventions, the free coordinate path integral is 

(y r 



Vy exp 



(4vrT) " 



(19) 



/o V 4 , 

All other free path integrals are normalized to unity. 
Scalar Field Theory with a self-interaction potential U((p) 

For a (real) scalar field with a self-interaction U(<p), the path integral represen- 
tation of the one-loop effective action analogous to eq.( 0) reads 



I [°°dT * T 
2Jo T 



J Vx exp 



T x 2 



+ U"{cj>{x)) 



(20) 



For example, for U ((f)) = ^(f) 3 one has an interaction exponential 



exp 



A / dT<j>{x{r)) 



leading to a scalar vertex operator such as in eq.( ||]). Wick-contraction of 
N such operators, using 



{y >1 (n)y»(T 2 )) = -g^G B (T 1 ,T 2 ) 

then leads to eq. ( [^) . 

Photon Scattering in Quantum Electrodynamics 



(21) 



We have already discussed photon scattering in scalar QED. Eq.( ||) generalizes 
to the Dirac fermion loop as follows (see e.g. H): 



T[A) 



dT 



T>x Dtp exp 



dr -x 1 + -^+ieAaX tM -ie^F^ 
\4 I 



(22) 

The "bosonic" coordinate path integral T>y is identical with the scalar QED one. 
The additional "fermionic" Dip - integration is over anticommuting Grassmann 
- valued functions, obeying antiperiodic boundary conditions, ip^(T) = —ip^(0). 
We may think of this double path integral as breaking up a Dirac spinor into 
a "convective part" and a "spin part". The global factor of (—2) compared to 
the complex scalar loop case is due to the different statistics and degrees of 
freedom. The photon vertex operator acquires an additional Grassmann piece, 



V A 



dr (e^ + 2ie^k v iP v Y k ' x{T) 



(23) 



The Grassmann path integral is performed using the Wick contraction rule 
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(f(ri)f(r 2 )) = \g» u G f {t u t 2 ), (24) 

G f {t 1 ,t 2 ) = sign(ri,r 2 ). (25) 

Before Wick-contracting two Grassmann fields, they must be made adjacent 
using the anticommutation property. The worldline action has a "worldline 
supersymmetry" , namely 

S v x^ = -2r/^, 8^ = tjx 11 , (26) 

with some Grassmann constant rj. It makes therefore sense to combine y and 
ip into a superfield. Defining 

X» = x» + V29^ = x%+Y», D = -^-6^-, Jd69 = l, (27) 
eq.( p2|) can be rewritten as a super path integral, 



!.'■[.. 1] 2 I ^e" m2T / VXexpl- jJlT Jd6 



00 dT 

-t J { 



--X D 3 X-ieDX fl A lx (X) 



(28) 

The two Wick contraction rules are then combined into a single rule for the 
field Y, 

(Y^T 1 ,6 1 )Y^T 2 ,e 2 )) = -g^G(r 1 ,6 1] r 2 ,e 2 ), (29) 
G(n,9i;T2,9 2 ) = G B ir l ,T 2 )+e x e^G F {jx^). (30) 
The photon vertex operator eq.( |2^) for the spinor loop is rewritten as 

V A = - [ T dTd9e.DX^exp[ikX]. (31) 
Jo 

This superfield reformulation makes no difference for the final integral repre- 
sentations obtained, but is often useful for keeping intermediate expressions 
compact. 

Photon Scattering in a Constant External Field F 

An additional background field A^{x), with constant field strength tensor F^, 
can be completely taken into account by changing the worldline Green's func- 
tions, and the path integral determinant [3^, 36, 37]. The Green's functions 
become 



Gf{ti,t 2 ) = G F 



12- 



Q—ieFTG bi2 

cos(eFT) 



(32) 
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These expressions should be understood as power series in the field strength 
matrix F (note that eqs.( |32|) do not assume invertibility of F). The background 
field breaks the Lorentz invariance, so that the generalized Green's functions 
are nontrivial Lorentz matrices in general. Accordingly, the Wick contraction 
rules eqs,(^l|),(p4|) have to be rewritten as 



(r(nW(T2)) = \Q^{n,T 2 ). (33) 

Again momentum conservation leads to the freedom to subtract from Qb its 
constant coincidence limit, 

Gb(t,t) = -L^(eFcot(eFT) - ~) (34) 
To correctly obtain this and other coincidence limits, one has to set 

sign(r, r) = 0, 

and adopt the general rule that coincidence limits must always be taken after 
derivatives. It follows that, for instance, 

G b (t,t) = 0, G 2 b (t,t) = 1. (35) 

The change of the path integral determinant eq.( |l^) induced by the external 
field is 



(4ttT)' 



(4vrT)" 



D_ 

' 2 



(47tT)" 

(4vrry 



D 

"~det" 



_D 

2 det" 



sin(eFr) ' 

eFT 
tan(eFT) 



eFT 



(Scalar QED) (36) 
(Spinor QED) (37) 



Scalar or Spinor Loop Contribution to Gluon Scattering 

The path integrals for the external gluon case differ from eqs.( |9|),( ^2|) only 
by path-ordering of the exponentials, and the addition of a global colour trace. 
For a given iV-gluon amplitude, this trace factors out as tr(T ai • • • T ajv ), where 
the T ai are the colour matrices carried by the gluon vertex operators, eq.( ||). 
The path-ordering leads to ordered r - integrals / Ilila 1 ^ r i^( T i — T i+i) such as 
in eq.( ||) (translation invariance in r is used for setting = 0). 
In the fermion loop case, the worldline Lagrangian eq.( ^) now contains a 
term ^[A^, A u ]ip u which, in the component formalism, forces one to introduce 
an additional two - gluon vertex operator besides eq.( |23|) [23]. This is not 
necessary in the superfield formalism, where the super vertex operator eq.( ^) 
remains sufficient. Here the only change is that a suitable supersymmetric 
generalization of the above 9 - functions is needed. This is 
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0(f y ) = 0(n - Tj) + 0i6j5{Ti - Tj), 



(38) 



where fa = n — Tj + OiOj. The nonabelian commutator terms above are then 



generated by the 8 - function terms [52|. 



Gluon Loop Contribution to Gluon Scattering 

The treatment of the gluon loop case in this formalism is considerably more 
involved than the scalar and spinor loop cases. This is due to the fact that the 
consistent coupling of a spin-1 path integral to a spin-1 background requires 
the introduction of auxiliary degrees of freedom, whose contributions have to 
be removed later on. I will just write down the appropriate path integral, and 



refer the reader to [23, 36] both for its derivation and evaluation. It reads 



1 f°° dT 

= » lim / -^exp 

l C'^oo Jq 1 
xtrPexpj— J dr 



- 1) 



Vx- 
P 2 



-x 2 + igA^ + ^[(dr 



A 

C)8n V 



VipVip 



2igF ftv ip 



(39) 



The Grassmann path integral now appears both with antiperiodic ( "A" ) and pe- 
riodic ("P") boundary conditions. V denotes the path-ordering operator. This 
path integral actually describes a whole multiplet of p - forms, p = 0, . . . , D 
circulating in the loop; the role of the limit C — > oo is to suppress all contribu- 
tions from p > 2, and the contributions from the zero form cancel out in the 
combination f^ — fp- 

Pseudoscalar Yukawa Coupling 

We consider now a pseudoscalar <f> interacting with a Dirac fermion tp via the 
A'0075?/> ~ vertex. In contrast to the gauge path integral, the following path 
integral representation for the one-loop effective action induced for the pseu- 
doscalar field by a spinor loop was found only recently 



S Yps [x,ip,ip5] 



POO 


dT 


7° 


~f 




fx 2 


f * 









Vx Vip Vip 5 e"° Yps 



S Y ps [x,1p,1p5] 



,2 i \2 i2t 



(40) 



Besides y and ijj there is now a second fermionic field "05- It is a Lorentz scalar, 
and has the same correlator as ip: 



(il>5 (n) ^5(72)) = -G f {ti,t 2 ) 



(41) 
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Scalar Yukawa Coupling 



The case of the scalar - fermion coupling Xift^ifj differs from the pseudoscalar 
case only by the following change of the action, 



Sy P s — > Sy s = Sy ps + 2Xm4>(x) 



(42) 



For the combination of both couplings see , for the nonabelian case |29| , |30| . 
To conclude this section, let us mention that the representation of axial vector 
couplings on the worldline has been studied in [31, [32], 33]. A treatment of 



external fermions along the present lines is still lacking. 



1.3) Multiloop Generalization 

If one wishes to generalize Strassler's approach beyond one loop, obviously one 
needs to know how the one- loop Green's function Gb, defined on the circle, 
generalizes to higher order graphs. As a first step in the construction of such 
"multiloop worldline Green's functions" J38| | , we will ask the following simple 
question: How does the Green's function Gb(t\,T2) between two fixed points 
T\,T2 on the circle change, if we insert, between two other points r a and r^, a 
(scalar) propagator of fixed proper-time length T (fig. 7)? 



T T \ 




Figure 7: Change of the one-loop Green's function by a propagator insertion 

To answer this question, let us start with the worldline-path integral represen- 
tation for the one-loop two-point - amplitude in (ft 3 - theory, and sew together 
the two external legs. The result is, of course, the vacuum path integral with a 
propagator insertion: 



Vx / dr a / dn {(ft(x(T a ))(ft(x(T b ))} exp 
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T ™2 



r (2) _ / -m 2 T 

vac J T e 

(43) 

Here (4>(x(r a ))(j)(x(Tb))} is the x-space scalar propagator in D dimensions, 
which, if we specialize to the massless case for a moment, would read 



<<Kx(t q ))<K*(t 6 ))> 



r( 



i) 



Air 2 



(x a - x b y 



(f-D 



(44) 
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Clearly this form of the propagator is not suitable for calculations in our aux- 
iliary one-dimensional field theory. We will therefore again make use of the 
Schwinger proper-time representation, 



roc _ _ „ _ _ D 
(cP(x(r a ))cP(x(r b ))) = / dfe~™ T (4vrT)--exp 
jo 



(x(r a ) - x(r h )Y 
4T 



(45) 



where the mass was also reinstated. We have then the following path integral 
representation for the two - loop vacuum amplitude: 



r (2) 



o 1 Jo Jo Jo 



dr h 



x / Vx exp 



T dr i 2 (x(r a ) - x(r b )Y 
o T 4 4T 



(46) 



The propagator insertion has, for fixed parameters T,T a ,r b , just produced an 
additional contribution to the original free worldline action. Moreover, this term 
is quadratic in x, so we may hope to absorb it into the free worldline Green's 
function. For this purpose, it is useful to introduce an integral operator B ab 
with integral kernel 



B a b(n,T 2 ) = <5(n - T a ) - <5(ti - T b ) 5(r a - r 2 ) - S(r b - r 2 ) (47) 
{B ab acts trivially on Lorentz indices). We may then rewrite 

(x(r a ) - x(T b )) 2 = / dr\ / dT 2 x(T 1 )B ab (T 1 ,T 2 )x(T 2 ) . (48) 

JO JO 

Obviously, the presence of the additional term corresponds to changing the 
defining equation for G B , eq.( 18), to 



4 1) (n,r 2 )=2(r 1 |(^-^y 



T2) 



(49) 



After eliminating the zero-mode as before, this modified equation can be solved 
simply by summing a geometric series: 



(- 



B ab 
f 



d^ 2 f d 



+1 



dr ) T ^dr > 



drJ 



B 



d^ 2 B 



— (— ) » , 

T \drJ T VdrV 
(50) 



leading to [38 | 



G ( b ) (t 1 ,t 2 ) = G b (ti,t 2 ) + 



1 [G B (n,T a ) - GB(n,T b )}[G B (Ta,T 2 ) - G B (T b ,T 2 )\ 



T + G B (r a ,n) 



(51) 
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The worldline Green's function between points t\ and ti is thus simply the 
one-loop Green's function plus one additional piece, which takes the effect of 
the insertion into account. Observe that this piece can still be written in terms 
of the various one-loop Green's functions Gsij- However it is not a function of 
7*1 — 7*2 any more, nor is its coincidence limit a constant. 

Knowledge of this Green's function is not quite enough for performing two- 
loop calculations. We also need to know how the path integral determinant is 
changed by the propagator insertion. Using the Indet = trln - formula, this 
can again be calculated easily, and yields 



J Vy exp 


[ fT^y 2 (y(Ta)-y(T b )) 2 - 


Def (£ _ A*) » 


J 


Vy exp 




Det'(^r ? 



G 



Bab 



T 



(52) 

As usual, the prime denotes the omission of the zero mode from a determinant. 
To summarize, the insertion of a scalar propagator into a scalar loop can, for 
fixed values of the proper-time parameters, be completely taken into account 
by changing the path integral normalization, and replacing Gb by Gg . The 
vertex operators remain unchanged. 

In the case of a photon insertion, instead of eq.( ^4|) one has, in Feynman gauge, 



=2T 



x(r a ) ■ x(r b ) 



D_ 

4ir 2 



£-1 

2 1 



(53) 



([x(r a ) -x(n)Y 

The denominator of the propagator is treated as above, i.e. exponentiated and 
absorbed into the worldline Green's function; the numerator x a ■ x^ remains, 
and participates in the Wick contractions. 

The whole procedure extends without difficulty to the case of m propagator 
insertions. The determinant factor eq.( p^) generalizes to 



(detA^ 



(54) 



and the two-loop Green's function to 



<4 m) (Ti,T 2 ) = G B (n,r 2 ) 



A 



(m) 
kl 



III 

+ o ^ [ G B(n,Ta h ) ~ G B (Tl,T bk ) 
k,l=l 

Here A^ is the symmetric m x m - matrix defined by 



G\B(r ai , 7*2) - Gb{j\, v T2) 
(55) 







c- 

T 


-1 






2- 




Tki 


= Tk5ki 




Cki 


= G B (T ak , 





(56) 
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and Ti, . . . , T m denote the proper-time lengths of the inserted propagators. 
Note that our formulas give the Green's function only between points on the 
loop. This suffices, for instance, to calculate single spinor-loop contributions 
to photon scattering in QED, i.e. sums of diagrams as shown in fig. 5, or 
to the analogous (pseudo-) scalar scattering diagrams in models with Yukawa 
couplings (the generalization to the case of several fermion loops interconnected 
by photon propagators requires no new concepts J|4| ) . 

In <f> s - theory or Yang-Mills theory there are, of course, also contributions to 
the amplitude with external legs on the propagator insertions. The extension 
of formula involving points located on the insertion was given in [38| for 



the two - loop case, and obtained for the general case by Roland and Sato [42|. 
This knowledge then is sufficient to write down worldline representations for 
all 3 graphs which have the topology of a loop with insertions. According to 
graph theory, for the first few orders of perturbation theory such a loop, or 
"Hamiltonian circuit" , can always be found - all trivalent graphs with less than 
34 vertices do have this property. </> 4 graphs are reduced to trivalent graphs in 
the usual way by introducing an auxiliary field, which can be again represented 
on the worldline H] . Even for 4> 5 an auxiliary field must be introduced if one 
wants to assure that all diagrams are generated with the appropriate statistical 
factors jw]]. A multiloop generalization for Yang- Mills theory still waits to be 
constructed (some preliminary steps have been taken in p6[). 



Roland and Sato also provided a link back to string theory by analyzing the 
infinite string tension limit of the Green's function G^f^ of the corresponding 
Riemann surface, and identifying G^ with the leading order term of Q^ s ^ 
in the \ - expansion: 

Gf (m \ Zl ,z 2 ) ^° IgM (71,75,)+ finite . (57) 

To be precise, denotes here not the Gg of eq.( |55|) but the equivalent 

Green's function G { ™\t 1 ,t 2 ) - {t x ,t x ) - jG { ™ } (t 2 ,t 2 ). 
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2. Second Lecture: Examples 

2.1) One-loop QED Vacuum Polarization 

For a warm-up, let us recalculate the one-loop vacuum polarization tensors in 
scalar and spinor quantum electrodynamics |23| . 

A ) Scalar QED 



According to eqs. the one-loop two-photon amplitude in scalar QED 

can be written as 



xx^(T 1 )e ikl - x{Tl) x u {T 2 )e ik2 - x{T2) e- ^ dT \ ±2 (58) 
Separating off the zero mode according to eqs.(]i3|),(|l7|), one obtains 



T^ihM = -(2n) D 5(k 1 + k 2 )e 2 ^~ mT drj dr 2 

JO -L Jo Jo 

x y , py^ t (n)e ifcl S/(Tl) ^(r 2 )e ifc2 y(T2) e"/o TdT ^ 2 . 



(59) 



We use the energy-momentum conservation factor (2n) D 5(ki + k 2 ) for setting 
ki = — k 2 = k, and then omit it. We now need to Wick-contract the two photon 
vertex operators, starting from the basic rule 

(^(ri)^(r 2 )) = -g^G B (n,T 2 ) . (60) 

The rules for Wick-contracting expressions involving both elementary fields and 
exponentials are the following: 

1. Contract fields with each other as usual, and fields with exponentials 
according to 

(^(n)e^ (T2) ) = WfaWteykve**™ (61) 
(the field disappears, the exponential stays in the game). 

2. Once all elementary fields have disappeared, contraction of the remaining 
exponentials yields a universal factor 



3 ifci-yi ,,, e ik N -y N \ =exp 



, N 



(62) 
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For the case at hand this produces two terms, 

(^(ri)e ife -^ ri )y !y (r 2 )e- iA; y ( T2 )) = ^ u G m2 - /WG| 12 Je" fc2Gsi2 . (63) 
Now one could just write out Gb and its derivatives, 

G b {ti,t 2 ) = sign(r 1 -r 2 )-2 (ri ~ T2) (64) 

G B (n, n ) = 25(n-T 2 )-^ (65) 

and calculate the parameter integrals. It is useful, though, to first perform a 
partial integration on the term involving Gbi2 in either t\ or t 2 . The integrand 
then turns into 

- ^r|G| 12 e- fc2Gsi2 . (66) 

Note that this makes the gauge invariance of the vacuum polarization manifest. 
We rescale to the unit circle, Tj = Tiii,i = 1,2, and use translation invariance 
in r to fix the zero to be at the location of the second vertex operator. We have 
then 

G B (n,Q) = Tu 1 (1-u 1 ),G b (t 1 ,t 2 ) = l-2«i. (67) 

Taking the free determinant factor eq.( |l9|) into account, and performing the 
global proper-time integration, one finds 



°°^e- m 2T (4vrT)-T T 2 /* du(l - 2u) 2 e~ Tu ^ k2 



(4vr)^ 



D 



i 2 Jo L 



--2 



(68) 



It is easy to verify that this agrees with the result reached by calculating the sum 
of the corresponding two field theory diagrams in dimensional regularization. 

B) Spinor QED 

For the fermion loop, the path integral for the two-photon amplitude becomes, 
in the component formalism, 



-2(-ief 



■dT 



-m 2 T 



Vx ( Vip ( dn [ dr 2 
J Jo Jo 



xixl + 2i$tyi ■ k^ kl ' xi (x v 2 + 2i^ip 2 • k 2 y ik2 - X2 e-fo dT (^+^ . 
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The calculation of T>x is identical with the scalar QED calculation. Only the 
calculation of T>i/j is new, and amounts to a single Wick contraction, 



-G 2 



F12 



^k 2 - k"k U 



t v k 2 - k"k v 



(70) 



Adding this to the bosonic result shows that, up to the global normalization, 
the parameter integral for the spinor loop is obtained from the one for the scalar 
loop simply by substituting, in eq.( |66|), 



r 2 

^B12 



G2 f~jl 
B12 ~ ^12 



(71) 



The complete change thus amounts to supplying eq.( ^) with a global factor 
of —2, and replacing (1 — 2u) 2 by — Au(l — u). This leads to 



(4tt)^ 



k^-g^k 2 ]v(2-—) f duu{l-u)\m 2 + u(l - u)k 2 
J 2 Jo L 



(72) 

again in agreement with the result of the field theory calculation. 
The qualitative features of this calculation generalize to the ./V - photon am- 
plitude. It is obvious that the Wick contraction of N photon vertex operators 
will lead to an integrand of the form 



Pn {Gsij , Gsij , GFij ) exp 



k<l 



Gskikk ■ h 



(73) 



with some polynomial Pn- What is remarkable is that the substitution above 
can be promoted to the following general substitution rule. As Bern and 
Kosower have shown, it is always possible to remove all second derivatives Gsij 
appearing in the integrand by a suitable chain of partial integrations. Once this 
has been done, the integrand for the spinor loop case can be obtained from the 
one for the scalar loop by simultaneously replacing every closed cycle of Gb's 
appearing, say G Bixi^G Biiiz ' ' '^-Biniu by its "supersymmetrization" , 



G Bixi^G Bi 



G Bi r , 



G Bi\iiG Biiiz ' ' ' Gb 



Gpi x i 2 Gfi 



Gfi 



(74) 

(up to the global factor of —2). Note that an expression is considered a cy- 
cle already if it can be put into cycle form using the antisymmetry of Gb 
(e.g. GsabGBab = —GBabGBba)- This rule may be understood as a remnant of 
worldsheet super symmetry [p|. 

The close connection between the scalar and fermion loop calculations may 
appear surprising, as it appears to have no analogue in standard field theory. 
It arises because the treatment of fermions in the worldline formalism (and in 
string theory) does not correspond to the usual first order formalism, but to 
a second order formalism for fermions. The Feynman rules for spinor QED in 



the second order formalism (see [53 1 and references therein) are, up to statistics 
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and degrees of freedom, the ones for scalar QED with the addition of a third 
vertex (fig. 8) f\ Note that those rules also apply to the external fermion case. 



a 



= +ig ( (*Jt - *$) - zo ap (k A )p ) 




k^-m 2 



= -ig 2 rf p 



Figure 8: The second order Feynman rules for fermion QED. They correspond 
to the ones for scalar QED with a third vertex added. This vertex involves a^ u = 
| [7^, j u ] and corresponds to the ip^F^ v ip v - term in the worldline Lagrangian. 
For the details and for the nonabelian case see [53]. 



2.2) Photon Splitting in QED 

Next let us consider an application involving constant external fields, namely the 
one-loop three photon amplitude in a constant magnetic field. This amplitude 
gives rise to photon splitting, a process of potential astrophysical interest whose 
first exact calculation was performed by Adler in 1971 [54]. 
We will again do both the scalar and the spinor QED cases. This amplitude is 
finite, so that we can work in D = 4. First we have to specialize our formulas 
for the general constant background field to the pure magnetic field case. We 
choose the B - field along the z - axis, and introduce matrices I03 and I12 
projecting on the t, z - and x,y - planes, so that 

^1 thank A.G. Morgan for providing this figure. 
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We may then rewrite the determinant factors eqs.(^) as 



(75) 



det 2 



det 2 



sin(e-FT) 



eFT 
tan(eFT) 



(eBT) 



eFT 



The Green's functions eq.(|32|) specialize to 



sinh(eST) 

(eBT) 
tanh(e£T) 



(76) 
(77) 



Gb(ti,t 2 ) = Gbi 2 Io3 



T cosh(zG , si2) — cosh(z) 



+— 



2 zsinh(z) 
T ( sinh(zGsi2) 
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2z V sinh(z) 



G B 12 *F 



a , \ A T , sinh(zGBi 2 ) T /cosh(zG jB i 2 ) 1\ „ 

smn(zj \ smn(z) z. 



Gb(ti,t 2 ) 
Gf{ti,t 2 ) 



f, j .J* ^cosh(zG B i2) ^ T zsinh(zGsi2) 

<-TB12l03 + 2 0l2 „ . : , , Ul2 + 2 — . —. r — 



Tsinh(z) / rsinh(z 

^ T , ^ cosh(zG B i 2 ) T sinhjzGsia) . * 

cosn(,zj cosn(z) 



»F 



(78) 



We have now introduced z = eBT, and F = -g. In writing Gb we have al- 
ready subtracted its coincidence limit, which is indicated by the "bar". The 
coincidence limits for Gb and Gf are also needed, 



Gb(t,t) 
Gf(t,t) 



coth(z) - - ) E 



tanh(z)zF 



(79) 



To obtain the photon splitting amplitude, we will use these correlators for 
the Wick contraction of three vertex operators Vq and V12, representing the 
incoming and the two outgoing photons. 
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A ) Scalar QED 



The calculation is greatly simplified by the peculiar kinematics of this process. 
Energy-momentum conservation ko + k\ + k 2 = forces collinearity of all three 
four-momenta, so that, writing —ko = k = ton, 



k\ = —k, k 2 = —k; k 2 



kj 



k 2 



k ■ k\ = k ■ k2 = ki ■ k 2 = 0. 



(80) 



Moreover, Adler [54] was able to show on general grounds that there is only 
one non-vanishing polarization case. This is the case where the magnetic vector 
k x e'o of the incoming photon is parallel to the plane containing the external 
field and the direction of propagation k, and those of the outgoing ones are 
both perpendicular to this plane. An appropriate choice of £0,1,2 leads to the 
further vanishing relations 



£1,2 • £0 = £1,2 • k = £i j2 • F = 



(81) 



This leaves us with the following small number of nonvanishing Wick contrac- 
tions: 



2 rT 



ii 1 ■■ 

<VbViV2> = — *n / dnexp ~2J u)iLOjng B ijn b\Qb\%^i 2j^ £ o£BOi« (82) 

i=0 L i,j=0 J 8=0 



For compact notation we have defined Qq = u;, 0)1,2 = —Mis- 
performing the Lorentz contractions, and taking the determinant factor eq.( ^) 
into account, one obtains the following simple parameter integral for this am- 
plitude: 



C sca i[w, ui,U2,B] 



m" 



8UJUJ1UJ2 Jo 



-m 2 T 



dTT- 



i=0 



^WiCosh^Gsoi) ex P^~^X! G Bij + 



z 2 sinh {z) Jo 
2 



1 



dn dT2 Gbyi 

T cosh(zGBij) 



i,j=0 



2z sinh(z) 



(83) 



Translation invariance in r has been used to set the position tq of the incom- 
ing photon equal to T. For better comparison with the literature, we have 



normalized the amplitude as in [54]. 
B) Spinor QED 

For the fermion loop case, we will now make use of the superfield formalism. 
As in the case without a background field, Qb and Qf can be combined into a 
super propagator, 



£(t-i,0i;t 2 ,6>2) = Gb\2 + 0i6 2 Gf 



12 



(84) 
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This allows us to write the result of the Wick contraction for the spinor loop 
in complete analogy to the scalar loop result eq. (jS2|) , 



2 f T f ri 2 

{V Q V 1 V 2 )=q{ / dn ddi exp ujibjjnQijn \ E\F>\D 2 Q\ 2 e 2 ^^£q-Pq^o 



i=0 



\j=0 



i=0 



(85) 



Performing the 6 - integrations and Lorentz contractions, we obtain 



Cspin 



rrr 



-m z T 



dTT- 



4ujiOiu>2 Jo z 2 sinh(z) 

2 



x j Q dn dr 2 exp j-- J>^- [G Bij + - ^ 



cosh(z)G , Bi2+wiW2(cosh(z) — cosh(zGBi2" 



+ 



sinh(z)cosh(z) 
LULU1LU2GF12 



UJl- 



cosh(zG SO i; 



uj 2 - 



cosh(zG , B0 2; 



cosh(2:) 



sinh(z) £ sinh(z) 

sinh(zG , soi) (cosh(2:) — cosh (26^02) J — (1 



(86) 



A numerical analysis of this three-parameter integral has shown it to be in 
agreement with other known integral representations of this amplitude fl37[ . 
However, the method of calculation improves in several respects over standard 
field theory methods. In particular, it allowed us to make use of the vanishing 
relations for the kinematic invariants essentially on line one of the calculation 
to reduce the number of terms. 

2.3) The Two Loop QED (3 functions 

We proceed to the two-loop level, and to a recalculation of the two-loop (3 - 
function coefficients for scalar and spinor QED. In field theory, the fermion 
QED calculation requires consideration of the diagrams shown in fig. 9, while 
for scalar QED there are some more diagrams involving the seagull vertex. 
Of the possible counterdiagrams contributions those from electron wave func- 
tion and vertex renormalization cancel on account of the QED Ward identity, 
however mass renormalization must be taken into account. The worldline for- 
malism applies only to the calculation of the bare regularized amplitude; the 
renormalization has to be performed in field theory. 

In the following calculation, which is a variant of the one presented in [39|, 
it will not be necessary to distinguish between individual graphs. We will 
again couple the two-loop path integral to some constant field background, 
and write down the effective action induced for the background field. The (5 - 
function coefficient will then be extracted from the divergence of the Maxwell 
term ~ F^F^, calculated in dimensional regularization. Of course we are 
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Figure 9: Diagrams contributing to the two-loop QED vacuum polarization 

free to choose any background we wish; we will therefore impose the condition 
F^ V F V \ ~ <5^, which leads to some simplifications. 

A) Scalar QED 

Applying the formalism developed in the first lecture to this case, we obtain, 
for the scalar QED case, immediately the following parameter integral for the 
two-loop effective Lagrangian: 



,(2) 
""seal 



( 4*r»(4) 



xdet 2 



2 ' Jo i Jo Jo JO 

sin(eFT)" 



eFT 



det 2 



T + 



n {pab + Qba 



(Va ■ Vb) 



(87) 



Here T and T denote the scalar and photon proper-times, and T a ^ the endpoints 
of the photon insertion moving around the scalar loop. The first determinant 
factor is the same as before, and represents the change of the free path integral 
determinant due to the external field; the second one represents its change due 
to the photon insertion, and generalizes eq.( |52| ) to the external field case. The 
two-loop Green's function eq.( [5l]) generalizes to 



QB{ruT b )}[QB{Ta,T 2 ) ~ GB(r b ,T 2 )] 

f + \ (g ab + Q ba 



G B J {Ti,T 2 ) = Gb(ti,T 2 ) + 



We use this Green's function for Wick-contracting the "left over" numerator of 
the photon insertion, which gives 



(Va ■ Vb) 



tr 



GBab + 



1 {GBaa - G Bab) (G Bab ~ ^Bbb) 

2 f + \(p ab + g ba 



(89) 



Care must be taken again with coincidence limits, as the derivatives should not 
act on the variables r a , T b explicitly appearing in the two- loop Green's function; 
again the correct rule in calculating (y a y b ) is to first differentiate eq.( |H) with 
respect to ti, t 2 , and set t\ = r a , t 2 = r b afterwards. 
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Next note that with F chosen as above Qsab + QBba is a Lorentz scalar, so that 



det 2 



-(£,/. 



Gbc 



ab 



+ Qba 



D_ 

2 



(90) 



and the T - integration can be trivially performed. Note that the term contain- 
ing d a b (see eq.( [32|)) vanishes in dimensional regularization upon performance 
of the T - integral (it corresponds to a massless tadpole insertion in field the- 
ory). After the usual rescaling r a ^ = Tu a ^ one expands the integrand in a 
Taylor expansion in F, and picks out the coefficient of tr (F 2 ). This gives 



42l^]=tr(F 2 )-^ 
(4tt) 

1 -R 
J^GBab 2 + 



(4) r d ^ 



' ^ /-"'// -m 2 Trpi-D 

1 — — 
GBab 2 + 



2 J Jo T 
D-4 4 



2 3(D - 2) 3D 



l ,i 

du a / dltfc 
JO 

28 8 - 



3D D-2 



2-± 



GBab" ' 

(91) 



The scalar proper-time integral gives the usual global T(4 — D), and the re- 
maining parameter integral produces Euler Beta - functions: 



/ du a [ du b G Ba b n 2 = f du a [u a (l - u a )] n ° = J3(n+l-^,n + l- ^Y 
JO Jo Jo v 2 2 ' 

(92) 

Expanding the result in e = D — 4, the 4r ~~ terms cancel, as expected. The i 
- term becomes 

4S,[^] - ^77^^ • (93) 

So far this is a calculation of the bare regularized amplitude. It must still be 
supplemented with a contribution from one-loop mass renormalization, which 
involves only one - loop quantities (see |SS(l ), 

A^HW^to-I^^. (94) 
Adding up both contributions 

e (47t) 



and extracting the /3 - function coefficient in the usual way (see e.g. [p5|), one 
finds the known result [56|, 

/?s ( 2i(«) = ^ • (96) 
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B) Spinor QED 



For the spinor loop in the superfield formalism one obtains again parameter 



integrals formally analogous to eqs.( |87|))( |89|) : 



dT 



1 i L> 

o T 1+ ~ 



xdet 2 



tan(eFT) 
eFT 



det 2 



-m 2 T 



1 



c?T / dr a dTb I d9 a d9b 
o Jo 



+ Q 



ba 



(-D a Y a ■ D b Y b ) 



(97) 



(-D a Y a ■ D b Y b ) = tr 



n n A 1 D a (QBaa — QBab)D b (QBab ~ Qsbb) 
D a L>byBab + 77 Z —7^ — 



(98) 



The further calculation also parallels the scalar case. The only point to be 
mentioned is that it is permissible and convenient to perform the T - integration 
before the 6 - integrals. The equivalent of eq.( ^) becomes 



£$ n [F] = tr(F 2 



(4vr) D Jo 



CO (]rp 

_ — m~_t jry 

T 



-T^Bab 2 + 




D 



4 + 



+ 



3L> 3(L> - 2) 



3D D-2 



Gbo 



This yields 



,(2) , , , 3a 



•^spin [-^1 



F. U F^ + O(e ) 



(4vr 



2 r v 



Mass renormalization contributes a [|39| 

5 „m r w, 4 a 2 



yielding a total of 



dm spmL 



,(2) 



e (47T 



F nv F^ 



2 r v 



1 



The (3 - function coefficient becomes 



(99) 



(100) 



(101) 



(102) 



a 

2^2 



(103) 



which is the classical Jost-Luttinger result [57]. 

I have chosen to show the most straightforward, though not the most efficient 
versions of this calculation. Both in the scalar and the spinor QED case the 
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2-loop f3 - function calculation can be further trivialized by suitable partial 
integrations in r a . The absence of a sub divergence for the Maxwell term then 
becomes manifest, regularization is only needed for the trivial global T - inte- 
gration, and all T a ,T b - dependence disappears from the integrand before the 
T a, T b ~ integration. This is discussed elsewhere |?9], Q. 

2.4) The Three — Loop Scalar Master Integral 

Finally, let us have a look at the simplest example of a three - loop parameter 
integral calculation in this formalism. This is the one where the integrand 
consists just of the bosonic three-loop determinant factor, eq.( |54]) with m = 2. 
In dimensional regularization, it reads 

rEP) = J o °°§e- m2T T e -l D I(D), (104) 



I{D) = / df 1 df 2 / dadbdcdd 
Jo Jo 



C 2 

(Ti + GBab){T2 + Gscd) — 



(105) 



y-. rp 

Here Xi^ = -m- denote the proper-time lengths of the two inserted propagators 
in units of T, and C = G Bac - G Bad - G Bbc + G Bbd . 

This is the most basic integral appearing if one applies this formalism to the 
calculation of three - loop renormalization group functions in abelian field the- 
ories, e.g. in QED |^3|] or the Yukawa Model. As far as the quenched QED j3 
- function is concerned, the general parameter integral appearing at the three - 
loop level still has a very similar structure: It has the same universal denomina- 
tor, possibly with a different power, and a numerator, which is again expressible 
in terms of the functions G Bab , G Bcd , C and their derivatives. 
In writing eq.( |104j ) we have already rescaled to the unit circle, and separated 
off the electron proper-time integral. This integral decouples, and just yields 
an overall factor of 

°° dT e -m*T T 6-lD = r( Q _ 3 ) m 3D-12 „ _£ 

o T K 2 1 3e 

The nontrivial integrations are Jq 1 da db dc dd = J abcd , representing the four prop- 
agator end points moving around the loop (fig. 10). 

This fourfold integral decomposes into 24 ordered sectors, of which 16 constitute 
the planar (P) (fig. 10a) and 8 the nonplanar (NP) sector (fig. 10b). Due to the 
symmetry properties of the integrand, all sectors of the same topology give an 
equal contribution. The integrand has a trivial invariance under the operator 
~§a + M + Jc + 33' wri ich just shifts the location of the zero on the loop. 
As a first step in the calculation of 1(D), it is useful to add and subtract the 
same integral with C = and rewrite 



I(D) = I sing (D) + I rcg (D), (106) 
I S m S (D)=J dfx df 2 J [(Ti + G Bab )(f 2 + G Bcd )\ 2 (107) 
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Figure 10: 10a 



10b 



-fsing(-D) factorizes into two identical three-parameter integrals, which are ele- 
mentary: 

r r°° r 1 r i~2.} 2 [2R(2 — — 2 — —)l 2 

I S ^(D) = {f o dT J o du[T + u(l-u)} 2 } =[ [ \ 2 _ 2 2) 

The point of this split is that the remainder I reg (D) is finite. To see this, set 
D = 4, expand the original integrand in — ^ — , and note that for all terms 
but the first one the zeroes of Gsab (Gscd) at a ~ 6 (c ~ d) are offset by 
zeroes of C 2 . Only the ^ - pole is required in renormalization group function 
calculations, so that one can set D = 4 for the calculation of I reg (D). The 
integrations over Xi,T2 are then elementary, and we are left with 



-^reg (4) 



4 / 



c 2 



1 
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AG Bab G\Bcd ' G,Ba£> Gficd 



(108) 



For the calculation of this integral, observe the following simple behaviour of 
the function C under the operation D ao = 



a I a 

dr a ' dr b " 



D ab C = ±2xnp, D 2 ab C = 2^8 ac - Sad - 5bc + &bd 



(109) 



where xnp denotes the characteristic function of the nonplanar sector. From 
these identities and the symmetry properties one can easily derive the following 
projection identities, which effectively integrate out the variable C: 



J /(C) GBab, 



G 



Bed) 



NP 



- A da 
Jo Jo 

f(C,G Ba b,G Bcd ) = -4 / da f dc 

Jo Jo Jo 



dc(a — c)/^— 2c(l — a), a — a 2 , c — c 2 ^j , 

2c(l-a) 



dc/(c, 



2 2 

a — a , c — c 



(110) 



Here / is an arbitrary function in the variables G,GBab,GBcd, and J C dCf 
denotes the integral of this function in the variable C, with the other variables 
fixed. The integrals on the left hand side are restricted to the sectors indicated. 
For / the integrand of our formula eq. (|108]), we have 
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c 



dCf 



C 



+ 



GBabGBcd C 

4 



ln( 



C 2 



+ 



VGBabGBcd 



^GBabGBcd 

arota„h(I 
V2 



C 



2 VGBabGBcd 



(111) 



Inserted in the second equation of ( |110| ) this leaves us with three two-parameter 
integrals, of which the first one is elementary. Applying the substitution 



y 



c(l - a) 
a(l-c) 



to the second integral, and 



2 c (! - a ) 



y 



a(l-c) 

to the third integral, those are transformed into known standard integrals, tab- 
ulated for instance in [58]. The result is 

/ / = 12C(3) - 8C(2) . 
Jnp 

The calculation in the planar sector is elementary, and I just give the result, 

^/ = 4C(2)-4 . 
Putting the pieces together, we have, up to terms of order O(e ), 



rE(£>) 



m 



3D-12 



r(6 - -d) 



2B(2 



D 
2 ' 



2-§)i2 



D 



12C(3) - 4C(2) - 4 
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